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Abstract 

An ecological flow network is a weighted directed graph in which nodes are 
species, edges are "who eats whom" relationships and weights are rates of 
energy or nutrients transfer between species. AUometric scaling is a ubiq- 
uitous feature for flow systems like river basins, vascular networks and food 
webs. By "ecological network analysis" method, we can reveal the hidden 
allometry directly on the original flow networks without cutting edges. On 
the other hand, dissipation law, which is another signiflcant scaling relation- 
ship between the energy dissipation (respiration) and the throughflow of any 
species is also discovered on the collected flow networks. Interestingly, the 
exponents of allometric law [t]) and the dissipation law (7) have a strong con- 
nection for both empirical and simulated flow networks. The dissipation law 
exponent 7 rather than the topology of the network is the most important 
ingredient to the allometric exponent 77. By reinterpreting t] as the inequal- 
ity of species impacts (direct and indirect influences) to the whole network 
along all energy flow pathways but not the energy transportation efficiency, 
we found that as 7 increases, the relative energy loss of large nodes (with high 
throughflow) increases, 77 decreases, and the inequality of the whole flow net- 
work as well as the relative importance of large species decreases. Therefore, 
flow structure and thermodynamic constraint are connected. 
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1. Introduction 



Ecosystem is a thermodynamic system driven by energy flows which origin 
from the sunhght a nd will be consumed by living organ i sms a nd dissipated 
to the environment (jOdunu . 1 19881 . Il983l : IStrakraba et all 119991 ) . In this sys- 
tem, hundreds of species interact each other and connected together by prey- 
predator interac tions to forrn . an entangled complex network which is alway s 



called food web( lPimml . 120021 : ICohen et al.l . Il990l : iPascual and Dunnd . l2005l ) 



In the past dec ades, some remarkable common patterns had been found in bi 
nary food webs(|Wi 



1997: Krause et al. 



2000 



Cattin et al. 



liams et al.l . l2002l : ISugihara et al.l . 119891 : iBersier and Sugihara 



20031) and also been reprodu ced by models (IWilliams and Martinezl . 
2004 lAUesina et aD . I2OO8I ) successfully. Food web. 



^ _ as 

the backbone of ecosystem, can transport energy flow from the environ- 
ment to every species as its unique function to differentiate from other net- 
works. In thi s view, food web posses s es som e special features, such as low 
trophic le vels (I Williams and Martinezl . l2000f ). energy bottlenecks and domi- 
nator tree ilAllesina and Bodinil . 120041 ). 

AUom etric scaling i s a remarkable universa l law found in variou s flow 



networksdKleiberj, Il932l : IWest et al.l . Il997l Il999l : iBanavar et al.l . 119991). The 



network embedded in d-dimensional space possesses a scaling relation C oc 
A^, where A is the metabolism or input flows from the source to the net- 
work, C is the total "mass" or the summation of all individual flow rates 
in the network and rj = {d + l)/d is the allometric exponent. River basins, 
mammalian blood vessels or plant vascular s ystems are the flow networks in 
d = 2 and 3 dimensi o nal sp a ces resp e ctivelvilRodriguez-Iturbe and Rinaldol . 



I997I : iBanavar et al.l . Il999l ). iDreyen (120011 ) further tested this assumption 



in one dimension a l spa ce by a water slot with evenly distributed sinks. 



Garlaschelli et al.l (120031 ) generalize the scaling relationship to the spanning 



trees of food webs. By cutting "weak" links of the original food web, they 
calculated Ai and Ci for all species on networks and found a power law re- 
lationshi p with exponent around 77 = 1.13 for almost al l the food webs they 
collec ted (lGarlaschelli et al.l . l2003l : ICamacho and Arenad . 120051 : iFrank and Murrelll . 



20051 ). 



However, the mentioned food web studies always neglected the energy 
flow informati on as the weight of links which is already available in the 
empirical data( lBrown and Gilloolyl . |2003|). Studies of ecological flow net- 
works concern ing both "who ea ts whom" binary relationship ar id "in what 
rate" problem (lUlanowiczl . l2004j ) in ecology have a long history (IFinnl . Il976 : 



2 



Szvrmer and Ulanowiczl. 119871 : iHigashil . Il986l : iBaird and Ulanowiczl . Il989l : iHigashi et al. 



1993; Patten. 1981 



1982 )■ By incorporating input-output analysis and Markov 



1986!), ecologists have developed a system atic ap- 



chain method (lHigashi . 
proac h called "ecological network analysis" in ecological flow networks ( jUlanowiczl . 



2004 ) . They designed a set of systematic i ndicators to describe the flow struc- 



tures and global state o f whole ecosystem (iFath and Patted. ll999l:lFath et al 



200ll: lFath_and Pattenl. Il998l : iHannonl . Il973l : iLevind . Il980l : iHannonl . Il986 : 



Pattenl . Il985l : lUlanowiczl . Il986l Il997h . 



AUometric scaling as an important pattern for all flow systems should be 
applied to ecological flow networks. However, due to the limitation of the 
existing approach developed by Garlaschelli so far, the extension of allometric 
scaling laws to any flow network is still lack due to three major reasons: 1. 
The original n aethod can be only applied t o trees so that many edges must be 
cut (however, lAUesina and Bodinil ( 120051 ) extended this method to directed 
acyclic graphs); 2. Information on flows and weights are never considered in 
previous works; 3. The ecological meanin g of allometric exponen t rj should 
be re-considered for general flow networks (IZhang and Gud . l2010l ). 

This paper applies the energy flow analysis method to calculate the key 
indicators of Ai and Q for all species on flow networks (Section [2]). The 
allometric scaling laws for 19 empirical ecological flow networks are shown 
in Section 13.21 Furthermore, we found that another exponent of the scal- 
ing relationship called the dissipation law (Section 13. 3p in this paper is the 
major ingredient to influence the allometric exponent. We have tested this 
hypothesis by a large number of numerical experiments both on empirical 
and simulated flow networks (Section 13. 4p . Finally, we re- interpret the ex- 
ponent in Section m as the indicator of inequality of species impacts, i.e., the 
concentration degree of the species impacts to the whole network on large 
species. 



2. Methods 

2.1. Review of Garlaschelli's Method 

To clarify the contribution of our method and its connection with the 
existing method, we review Garlaschelli's method for a hypothetic food web 
at flrst. 

Figure [T](a),(c) shows how Garlaschelli's approach can be applied to a 
hypothetical flow network (a) to calculate and Cj for each node. At 
flrst, a spanning tree (Figured] (c)) is constructed from the original network 
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(a) 




50 





(c) 
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(d) 



Figure 1: Comparison of different methods calculating the allometric scaling of a hypo- 
thetical ecological flow network. 

(a) is a hypothetical network (The letter in each node is its index) . The black node 

is the root, the numbers besides edges stand for flows, the dashed lines represent 
dissipations; (b) shows the Ai,Ci values being denoted inside and beside node i 

rpsr»pr*ti'^7pK7 hi\T nnr mptlinrl* { o\ is pt S"nannincf trpp r»f ( a\ tlip TinmKprs inftirlp anrl 



(Figured] (a)) by cutting edges. That way, each sub-tree rooted from any i 
can be viewed as a sub-system of the spanning tree. Ai is the total number 
of nodes involved in this sub-tree and Ci is the summation of AiS for each 
node in this sub-tree. Finally, the universal allometric scaling relationship 
of AiS and Qs, with an exponent aro und 1.13, was found for all food webs 
according to iGarlaschelli et al.l (|2003[ ). 



Garlaschelli's method was inspired by iBanavar et al.l ( ll999l )'s model to 
explain the Kleiber's law. The spanning tree is simply Banavar's optimal 
transportation network. Thus, energy flows into the whole system from the 
root along links of the network to all nodes. Suppose that each node would 
consume 1 unit of energy in each time step. A flux with 1 unit representing 
the energy consumption by each node should then be added to the orig- 
inal spanning tree Figure Wic). In Figure HTd), the ener gy dissipation by 



each node is added as a dotted line (IBanavar et al.l . Il999f ). As a result, A 



of each node is the throughfiow of this node. Ci is the total throughfiows 
of the sub-tree rooted from i. Essentially, calculation of allometric scal- 
ings u sing Garlaschelli ' s app roach is based on this weighted flow network 
mode hlZhang and Guol |2010| ). 



2.2. Ecological Network Method 

We will extend Garlaschelli's method directly on the original weighted 
network without cutting edges (e.g. Figure [T]). But the key question is how 
to calculate A^ and Ci for general flow network? 

According to the flow network interpretation of Garlaschelli's method in 
Figure [1] (d), Ai is the energy flux intake by node i which is equals to the 
total out flows from i due to the flow balance condition. Thus, this concept 
can be extended to any flow network by defining Ai as the throughfiow of 
node i. However, defining Ci is not as simple as Ai because we don't have 
the sub-system concepts if the considered network is not a tree. 

To understand what does Ci mean in Garlaschelli's method based on the 
flow network picture (Figure [T](d)), we suppose to do the following hypothetic 
experiment. Assume a large number of particles are flowing along the network 
(Figured] (d)), and all particles whenever passing by any node, say b, will be 
attached a label, say "b". This trace marker will not be erased forever unless 
the particle flows out of the network. Then, we found Cb in Figure [1] (c) is 
just the total number of labeled particles by "b" still being trapped in the 
whole network. This trace marker experiments can be also applied to other 
nodes independently and separately by attaching different labels so that all 
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CjS can be calculated by counting the number of particles who used to pass 
node i. 

Actually, this understan ding can be extended to any flow network no 
matter if it is a tree or not flZhang and Guol . l20ld ). Although counting the 
number of labeled particles in the real network is difficult, we can do this 
calculation directly by Markov chain technique thanks to the ecological net- 
work analysis method developed by Patten et al. As long as the flow network 
is in the steady state so that all flows distributed on edges are stable, a flxed 
Ci value can be calculated according to the flow structure. 

Suppose the flux matrix of the original network is F, in which each entry 
fij stands for the flux from i to j. Two special nodes and + 1 (A^ is the 
total number of species) representing source and sink are contained in this 
matrix as the flrst(last) in columns (rows). We deflne Ai of i as. 



N+l 

A = 5^/ii,VzG[l,Ar] 



(1) 



That is the throughflow of species i in ecological network analysis ( jUlanowiczl . 
1997h . 

To calculate Ci, we should convert the original flux matrix F into a 
Markov chain in which each element is deflned as rriij = fij/{Yl!k=i fik) for all 



1 < i < N . Notice that the Markov chain is normalized (i.e. ^ij = 1) 

only if the original flux matrix is balanced, i.e., 

N+l N+l 

E/^^ = E/^^'^^^[i'^]- (2) 

i=0 j=0 



This requirement is always satisfled by most empirical ecological net- 
works. The webs not satisfying this condition will be balanced by the method 



mentioned in Appendix Appendix A Thus, Ci can be calculated as 



N N 



C. = J:Y. for '''''' 



k=l j=l 



(3) 



where, Ujj is the elemen t in matrix f/, which is called fundamental matrix( lFath and Fatten! . 
19991 : lUlanowiczl . 120041 ) and deflned as. 
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U = {I-M)-\ 



(4) 



where, / is the identity matrix. According to the ecological network analysis 
method, Cj is just the total number of particles that used to pass node i (see 



Appendix Appendix B). 



Then, we can calculate A,-, C, values for all nodes in the flow network to 



test the following allometric scaling law, 

a oc Ai 



(5) 



Where, t] is the allometric exponent that will be mainly discussed in the 
following sections. 



3. Results 

3.1. Description of Data Set 

The 19 ecological flow networks in different habitats are studied(Table 
[1]). The networks are obtained from the online data basjll. Most of these 



networks in this database are from the published papers (lAlmunia et al.l . 11999 

ic3, 



Baird and Ulanowiczl. Il989l: iBaird etaP. Il998l: iHagvl. I2OO2I: luianowic: 



Monaco and Ulanowicz . 199?! : Christian and Luzkovich . 19991 ). In Table [H 



1986 



we list the name, the number of nodes [N) and the number of edges {E) 
of these networks. In which, the nodes are living species and also non- 
living compartments (e.g., DOC, POC), the weighted links are energy flows 
whose values vary in a large range because the units and time scales of the 
measurements are very different. The source node and sink node N + 1 are 
the "input" node and the combination of "respiration" and "output" nodes 
in the original data respectively. The dissipative flux of each node i is just 
the flow from i to N + 1 which can be read from the data directly. Most 
ecological flow networks are balanced already (condition Equation [2]), few 
imbalanced networks are balanced by the approach mentioned in Appendix 



Appendix A 



3.2. Allometric Scaling Law 

We found all the ecological flow networks possess the allometric scaling 
pattern signiflcantly as an example shown in Figure |21 Their allometric 



http://vlado.fmf.uni-lj.si/pub/networks/data/bio/foodweb/foodweb.htm 
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Table 1: Empirical Ecological Flow Networks and Scaling Exponents 
N and E are the total numbers of nodes and edges respectively, all networks are 

sorted in decreasing order of N 



bood web 


Abbre. 


AT" 

iV 


IT' 




V 


Rl 




7 




Florida Bay, Dry Season 


Baydry 


126 


2102 


1 


.010 


0, 


.995 


0. 


,915 


0. 


,949 


Florida Bay, Wet Season 


Baywet 


126 


2071 


1 


.020 


0, 


.994 


0. 


,917 


0. 


,953 


Mangrove Estuary, 
























Dry Season 


Mangdry 


95 


1462 


1, 


.010 


0, 


.997 


0. 


,978 


0. 


,983 


Everglades Graminoids, 
























Wet Season 


Gramdry 


67 


863 


1 


.030 


0, 


.999 


0. 


,973 


0. 


,997 


Everglades Graminoids, 
























Wet Season 


Gramwct 


67 


863 


1 


.020 


0. 


,999 


0. 


,977 


0. 


,998 


Cypress, Dry Season 


CypDry 


69 


639 





.998 


0. 


,996 


0. 


,957 


0. 


,949 


Cypress,Dry Season 


CypDry 


70 


554 





.998 


0. 


,996 


0. 


,967 


0. 


,949 


Cypress,Wet Season 


CypWet 


69 


630 





.997 


0. 


,997 


0. 


,965 


0. 


,988 


Mondego Estuary 
























-Zostrea site 


Mondego 


44 


401 


1 


.010 


0, 


,999 


0. 


,979 


0. 


,997 


St. Marks River (Florida) 


StMarks 


52 


349 


1 


.020 


0, 


,980 


0. 


,985 


0. 


,950 


Lake Michigan 


Michigan 


37 


210 


1 


.010 


0. 


,999 


0. 


,995 


0. 


,999 


Narragansett Bay 


Narragan 


33 


194 


1 


.010 


0, 


.991 


0. 


,813 


0. 


,942 


Upper Chesapeake 
























Bay in Summer 


ChesUp 


35 


203 


1 


.050 


0, 


.997 


0, 


.952 


0. 


,991 


Middle Chesapeake 
























Bay in Summer 


ChesMiddle 


35 


195 


1 


.040 


0. 


,996 


0. 


,851 


0. 


,761 


Chesapeake Bay 
























Mcsohalinc Net 


Chesapeake 


37 


160 





.994 


0, 


.997 


0, 


.985 


0. 


,985 


Lower Chesapeake 
























Bay in Summer 


ChesLower 


35 


163 


1 


.050 


0. 


,997 


0. 


,926 


0. 


,971 


Crystal River Creek 
























(Control) 


CrystalC 


22 


107 


1 


.040 


0, 


.997 


0, 


.959 


0. 


,995 


Crystal River Creek 
























(Delta Temp) 


CrystalD 


22 


83 


1 


.040 


0, 


.998 


0, 


.963 


0. 


,996 


Charca de 
























Maspalomas 


Maspalomas 


22 


82 





.956 


0, 


.966 


1 


.150 


0, 


.737 


Rhode River Watershed 
























- Water Budget 


Rhode 


18 


54 





.828 


0, 


.866 


1, 


.200 


0, 


.963 
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exponents r^s with the values of i?^ are listed in Table [U We see that the 
values of are larger than 0.9 for all food webs except CrystalC, CrystalD 
and Rhode whose scales are very small [N < 23). All of exponents rj fall into 
the interval [0.83, 1.05], most of them are larger than 1 a little. 

To test if the allometric scaling pattern is significant compared to random 
flow networks, we built a null model in which the numbers of nodes and edges 
are kept, all links are re-connected randomly and all flows on edges are also 
randomly assigned on the interval (0, fm] evenly, where fm is the maximum 
flux of the original network. From the inset of first row in Figure [21 we see the 
null model network doesn't show the allometric scaling law. We also try to 
compare the empirical flow networks with other null models. Some of them 
are obtained by keeping the topology unchanged but randomly assigning 
the weights, some of them are obtained by just shuffling the weights among 
edges, all the details of these null models are presented in the Appendix 



Appendix C 



Among these null models, we found the ones keeping weights information 
have similar allometric scaling exponent as the the empirical food webs. As 
a result, we know that it is the flow distribution but not the topological 
structure of the network that is the key ingredient to the allometric scaling 
exponents. However, to study all possible flow distributions which can affect 
the flnal allometric law is impossible because there are hundreds of flows 
which can be adjusted for an empirical food web. What is the most important 
aspect? We found the dissipative flux is the key to the allometric exponent. 

3.3. Dissipation Law 

In ecology, dissipation of a species has different forms such a s respiration, 



excre tion, egestion, natural and predatory mortality and so forth (IStrakraba et al 



19991 ). In our data, the dissipative flow is mainly respiration. We can un- 
derstand the dissipation of a species Di in an ecological flow network as the 
flows out of the network, i.e., Di = fi^N+i- It is comprehensible that this 
output flow increases with the total throughflow of the focus species. But it 
is not obvious that for most collected ecological flow networks the growth of 
dissipation along different species is slower than the growth of throughflow. 
A sub-linear relationship between dissipation and throughflow of each species 
is hold per se, 

A oc AJ, (6) 
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Allometric Scaling Law 




Figure 2: Allometric scaling law and dissipation law for Mondego ecological network and 
the null model (insets) 

There are only few points in the inset of the lower figure because many nodes in 
the null model are not balanced and their dissipations are set to zeros 
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Figure 3: The relationship between 7 and rj in original and adjusted ecological flow net- 
works 

The blue dashed line is the 7 — 77 curve for the original flow networks, and the solid 
red one is for the adjusted networks based on the same network structure and the 
dissipation law according to flow adjusted algorithm (see main text) 



where, the exponent 7 in Equation |6]is called dissipation law exponent which 
can be estimated from the data. Its values for different networks are also 
listed in Table [H We see all the 7 values are smaller than 1 except Maspalo- 
mas and Rhode networks. The lower plot in Figure |2] shows the dissipation 
scaling law for Mondego flow network as an example. 

3.4- Relationship Between 7 and rj 

Because both 7 and rj are indicators for the whole flow network, to see 
how these two numbers correlate each other, we can simply plot different 77s 
against 7s across all the collected empirical flow networks(see the blue dotted 
curve in Figure [3]) . Although a general trend that r} decreases with 7 can be 
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observed, it is not significant due to three major reasons. First. The number 
of sample points is too small to show clear relationship. Second, most of data 
points are concentrated in the circled area since all 7s and r^s are of similar 
values for all networks. Finally, the original exponents are noisy so that the 
rj values fluctuate in the circled area. 

To understand how rj depends on 7 for a given flow network structure 
as well as to avoid the problems mentioned in the previous paragraph, we 
invent a specific technique called "Flow Adjusting Algorithm (FAA)". This 
algorithm enables us to perturb the given flow network structure as little as 
possible and simultaneously observe how r] changes by tuning 7. Concretely, 
for a given original flow network F, we keep the network topology, the relative 
importance of each influx unchanged (That is, fji/J2jfji — f'jil Ylij f'jdi 
meanwhile adjust the flow distributions on each edge to obtain a new flow 
network F' such that: 

(1) The given dissipation law, i.e., D[ oc A'^ , is hold for every node i in 
F', where 7 is a given exponent which can be tuned. We will study how 7 
impacts the allometric scaling law; 

(2) The flux balance condition, i.e., Ylij f'ij — f'ji i^ust be kept for 
each node i; 

In this way, we can perturb the flow structure of the original network to 
obtain the expected dissipation law, and to observe how exponent 7 affects 
exponent rj. The details of the "flow adjusting algorithm" will be introduced 
in Appendix Appendix D Figure H] shows the dependence of r] and 7 on the 
perturbed networks by MAA based on the original Mondego flow network, 
randomized Mondego flow network (based on Mondego' s topology but assign 
flows randomly) and simulated networks by Niche mode KlWilliams and Martinezl . 



2OOOI ) . We observe the allometric scaling exponent decreases with the dis- 



sipation law exponent in a similar manner no matter the original network 
structures are. However, the concrete shapes of the curves between rj and 
7 change with the network structures. For networks generated by Niche 
model, rj decreases with 7 in a slower speed when the connectances are higher 
(blue triangles versus purple triangles and yellow diamonds versus green di- 
amonds). As a result, the dissipation law exponent but not the structure is 
the major feature to affect the allometric exponent. But we cannot conclude 
the topological structure has no influence on the allometric exponents, we 



will discuss this problem further in Appendix [Appendix E 



In Figure IH the red circle in the middle of the black curve stands for 
the combination of the dissipation law exponent 7 of the original Mondego 
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Figure 4: The relationship between 7 and rj by flow adjusting algorithm applying on 
Mondego ecological network, randomized Mondego flow network, and generated networks 
by Niche Model. 

The black and red curves are the original and randomized flow networks of 
Mondego respectively. The randomized flow network keeps the topological 
structure of Mondego network but assign flux randomly. Other dashed curves 
are for networks generated by Niche model with different number of nodes 
(100 or 50) and connectances (0 .4 or 0.1), as well as the random assigned 
flows! Williams and Martinez . 2000) • MAA is applied on all these flow networks to 
obtain the relationships between r] and 7. All the exponents are the average results 
of 50 random experiments. The red star shows the original position of dissipation 
law exponent (0.980) and allometric exponent (1.012) of Mondego network. While 
the red circle corresponds to the combination of the original dissipation law ex- 
ponent (0.980) and the adjusted allometric exponent (1.011) by FAA on Mondego 
network. 
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network F and the allometric exponent rj' of the adjusted flux matrix F' by 
the flow adjusting algorithm. While, the red star corresponds to the original 
exponents of Mondego network both for 7 and rj. We see these two points 
almost overlap together which means the original Mondego food web satisfies 
two conditions: 

(1) dissipation law is significant for the original 7; 

(2) the balance Equation [2] is obeyed exactly. 

However, for other empirical networks, the perturbed result of allometric 
exponent is not identical to the original one exactly on the given original 
dissipation law exponent because either the flux balance condition is violated 
or the dissipation scaling law is not significant (see Appendix Appendix F). 

By adjusting the fiows on all empirical fiow networks with the fixed orig- 
inal dissipation law exponent 7, we can eliminate noise in raw data because 
we conform the network to be fiow balanced and satisfy the given dissipation 
law. The red solid curve in Figure |3] shows the clear dependence of 7 and 77 
of all the empirical fiow networks. We see the decreasing trend of the blue 
dashed line standing for the original exponents is not as obvious as the red 
one representing the adjusted results. Therefore, the two conditions listed 
before cannot be satisfied by the some original ecological networks perfectly. 



4. Discussion 



4-1. Transportation Efficiency or Inequality of Species Impact? 

Previous studies explained the allometric scaling exponent 77 as the trans- 
portation efficiency of the network and is bound in between 1 (a star network, 
the most efficient tree) and 2 (a chain, the most inefficient tree). 

However, in our study, the allometric scaling exponent is not bound in the 
interval [1,2], therefore, we should give a new explanation for the exponent 
rj. The key problem is to understand the indicator Cj. 

In Banavar's model and Garlaschelli's method, they understand Ci as a 
cost of energy transportation for the sub-tree rooted from i. In this way, 
the energy fiows on the redundant links (loops or cross-leveled links) except 
the ones in the spanning tree are wasted. Nevertheless this interpretation 
can hardly generalize to fiow networks because (1) the wasted energy in the 
weighted networks can be measured as dissipation of each node but not the 
weight of edges; (2) all energy links should be considered because they all 
contribute to the whole fiow distribution. 
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According to the particle coloring experiment mentioned in Section [21 we 
can under stand Cj as the tot al impact of i to the whole network along all flow 



pathways( lVitali et al.l . 120111 ) because it is the total number of particles who 
used to pass i at least one time. Thus, as Cj increases, more other nodes will 
be affected by the particles used to pass z, the direct and indirect influences 
of i will increases. This understanding toward Ci can be extended to any 
flow network. 

As the species climbs up along the energy throughflow gradient Ai, its 
total impact Ci also increases with a relative speed r] according to the allo- 
metric scaling law Ci oc A^. Therefore, the important nodes (with larger Ai) 
may have much greater power (total impact Ci) in the networks with larger 
exponent rj than those networks with smaller exponents. For example, we 
have two networks with 4 nodes. They have the same throughflow distri- 
butions, e.g., Ai = {1,2,3,4}, but different exponents rji = 1 and ri2 = 2, 
therefore different d distributions Cf^ = {1, 2, 3, 4}, ^ = {1,4,9,16}. 
The most important node (with largest A^ = 4) in the second network may 
have much greater impact {C^^ = 16) to the whole network than the flrst one 
{C^^ = 4). Hence, the inequality of the species impact of the second network 
is larger than the flrst one. In short, the allometric exponent t] measures the 
inequality of species impact. 

This new interpretation is compatible with the previous one. For Gar- 
laschlli's spanning trees, the most inequable tree with a given root is a chain 
but the most equable one is a star. 

4-2. Dissipation and Inequality 

According to the new interpretation of allometric exponent and the dis- 
covery of the relationship between dissipation and allometry, we can obtain a 
whole picture: the networks may become more equable by dissipating more 
energy on larger nodes since the impacts of high flux nodes are weakened. In 
the networks with 7 > 1, the dissipation flow per throughflow of each node 
increases, the energy invested to the whole network decreases with the size of 
node, the network is more decentralized. On the other hand, if 7 is smaller 
than 1, the dissipating flux scales to the throughflow with a smaller relative 
speed. So the large nodes may input more energy on the whole network to 
obtain much more powerful impact on other nodes, the networks are more 
centralized. 

However, an interesting unexplained fact is the allometric exponents of 
empirical ecological networks are all close to 1. They are neither inequable 
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nor equable. We guess this exponent should be an optimal result by some 
unknown factors. This will be left for future studies. 

In summary, this paper generalize the universal allometric law to the 
ecological flow networks and discover that the major factor to the allometric 
exponent is 7, the dissipation law exponent. By reinterpreting allometric 
exponent as the inequality of species impacts, we build a connection between 
network structure and the thermodynamic constraint. This connection is 
very important deserving more attention. 
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Appendix. 

Appendix A. Balancing A Flow Network 

For most of empirical ecological flow networks, the flux matrix F is bal- 
anced which means J2f=o fji — "^^=1 /»i holds for each i G [1, iV]- However, 
some empirical networks and most artificial networks (e.g. random network) 
are imbalanced. Therefore, we should balance the given network F artificially 
so that the ecological network analysis methods can be applied. 

Suppose fji 7^ '^f=i fij node i. We can add an edge with the 

flux \w\, w = J2j^=o fji ~ '^j'=i fij connect node i to + 1 or 0. If w > 0, 
the direction of this artificial edge is from i to N + 1. Ifw<0, the direction 
is from to i. We can do this process for all nodes except and iV + 1 to 
balance the whole network. 

Appendix B. Explanation on Cj Calculation by an Example Net- 
work 

In this Section, we will explain why the number of particles labeled by 
node i can be calculated as Equation [3] through an example network (see 
Figure IB.Sp . 

Notice that the original network is imbalanced, we should balance it at 
first by the approach mentioned in the last section. Then, the balanced 
network can be converted to a Markov chain M simply through being nor- 
malized by the output fiow of each node as shown in Figure IB. 51 (c). We 
know any element rriij in Markov matrix M stands for the transfer proba- 
bility of a particle from i to j given that the particle locates i already and 
the whole network is in the steady state. And any element in matrix M* 
is the probability of a particle locating on i at first and transfer to j along 
all possible pathes after t steps. All these information is aggregated in the 
fundamental matrix U because, 

U = I + M + + ■ ■ ■ + M"^ = {I - M)-\ (B.l) 
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Figure B.5: An example network showing the calculations of Ai and Ci 
(a). The original flow network which is imbalanced; (b). The balanced network; 
(c). The corresponding Markov chain 

We should be careful to give an intuitive explanation on U because its 
element Uij does not stand for probability anymore although the element 
in each term M* is the transitional probability. We can write down the 
fundamental matrix for our example network (by ignoring the elements for 
source and sink N + 1), 



U^I + M + M^ + ---^ 



( 1 3/5 7/20 1/5 2/5 \ 

42/41 7/82 14/41 28/41 

12/41 42/41 4/41 8/41 

3/164 21/328 165/164 21/41 

V 3/82 21/164 1/82 42/41 / 



(B.2) 



Notice that the elements on the diagonal are larger than 1, so they cannot 
be interpreted as probability simply. 

Next, we will calculate the first passage flow to any node i. Here, Gj 
is defined as the number of particles passing i in the first time. If all the 
particles passed i arc labeled, then Gi is the number of particles passing % 
and unlabeled by it at each time. This quantity can be calculated as. 



N 



i=i 



Mi 



(B.3) 



For example, the first passage flow of node 4(^4) can be calculated as. 
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G4 = f: ^ = /o,i^ = 100 X -iZf- = 656/33 (B.4) 

^ M4,4 M4,4 165/164 

Actually, each term of Equation IB.3I is the first passage flow from node 
j to i, that is the number of particles that have visited j (in whatever time) 
and finally arrive at i in the first time along all possible flow pathways. By 
dividing by the term ua to derive first passage flow one av oids duplicate 



counting the particles who have visited i llHigashi et al.l . Il993l ). 

Because at each time, there are totally unlabeled particles will pass i 
and be unlabeled by "b", thereafter, Gi new "b" particles will be injected to 
the system and will flow to other nodes along pathways. Then, at given time 
step, there are totally Gi labeled particles just injected i, and ^^if GiM^ik} 
labeled particles injected one time step ago, and Yl^Zi ^i^fik} labeled par- 
ticles injected two time steps ago,...,^^~f^ GjMi-^i labeled particles injected 
t time steps ago, and so forth. 

Hence, the total number of labeled particles that flowing in the whole 
network at each time, being defined as Gi, is just the summation of the 
labeled particles being attached 1 time step ago, 2 time steps ago,..., and so 
on. Therefore, we can calculate Gi as: 



N 



a = ^ ^^(1 + M + + • ■ • + M°°; 



{ik} 



k=l 
N 



= J2GiU{ik} (B.5) 

k=l 
N N 

= ^ fojUji/Uii)Uik 

k=l j=l 

For example, G2 of node 2 in the example network is calculated as: 

5 5 
C2 = G2 ^ U2k = ((100 X 3/5 + 0)/(42/41)) «2fc = 125 (B.6) 

k=l k=l 



22 



Appendix C. Null Models 



To test if the allometric scaling law is a significant pattern in empirical 
ecological network, we designed four kinds of null models based on empirical 
networks. 

Null Model 1 (NMl): We only keep the total number of nodes and edges 
as the original empirical flow network, and build random connections, assign 
random weights for each edge. The weights are evenly distributed on the 
interval (0,/^], where fm is the maximum flux in the original network. In 
this way, the topology and flow distribution are destroyed. 

Null Model 2 (NM2): The connections are kept, the weights are randomly 
assigned for each edge. In this model, weights are also randomly sampled 
from the interval (0, fm]- In this way, only the flow distribution is destroyed. 

Null Model 3 (NM3): Keep the connections, shuffle the weights on edges. 
That is, we keep the topology and weights distribution but permute these 
weights on edges. In this way, the flow distribution is not changed but the 
correlations between flows are destroyed. 

Null Model 4 (NM4): Keep the weights, the number of edges, but ran- 
domly assign the weighted connections between any pair of nodes. In this 
way, the flow distribution is kept, but their correlations and the network 
topology are destroyed. 

For each original ecological flow network, we built four null models. The 
flux matrix F of the null models is imbalanced normally, then we should 
balance it by the approach mentioned in Section Appendix A after that we 
calculate their Ai and Ci for each node and derive the allometric scaling law 
pattern. 

Figure IC.6I shows the allometric scaling relationships for AiS and CjS for 
null models of Mondego network. From this flgure, we know NM3 and NM4 
have more similar pattern as the original networks than NMl and NM2, 
which means the flow distribution is more important than topology for al- 
lometric scaling. Although NM3 and NM4 have signiflcant scaling pattern, 
their exponents rjs are smaller than the one of the original Mondego net- 
work. Therefore, NM3 and NM4 cannot reproduce the main characters of 
the original network. 

Furthermore, we generate 50 networks for each null model on each em- 
pirical flow network collected. The average values of rjs are compared to the 
original ones in Figure [UTfl From this flgure, we can see that the exponents of 
NM3 and NM4 are more close to the original networks with less fluctuations. 
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Null Model 1 



Null Model 2 




Figure C.6: Allometric Scaling Patterns for Null Models based on Mondego Network 
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Figure C.7: Comparison for Allometric Scaling Exponents of Original Networks with Null 
Models 

All the empirical flow networks are sorted in the order of Table [T] along the hori- 
zontal axis(the left most network has largest number of nodes). The data points 
and error bars in all null models stand for the average numbers and standard 
deviations of exponents r] of 50 experiments. 



That means the weights information is more important than the structures 
and the weight correlation play a minor role on allometric exponents. Al- 
though NM3 and NM4 have similar exponents as the original networks, all 
their values are smaller than the ones of the original networks. 

Further studies on the dissipation laws of these null models can explain 
the patterns shown in Figure IC.7I From Figure IC.8[ we observe that most 
networks do not have obvious dissipation scaling law. However, we can ob- 
serve the obvious straight lines formed on the ceiling of data clouds for NM3 
and NM4, although lots of scatter points are below them. Actually, the arti- 
ficial balancing method can account for this phenomenon. Because in NM3 
and NM4, all the weights of links are not changed but the connections are 
destroyed so that the energy influx cannot balance with out flows for lots of 
nodes, the balanced flows (dissipations) are almost proportional to the origi- 
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Figure C.8: Dissipation Law for Null Models of Mondego Network 
The dissipation laws are calculated for the artificial balanced networks to compare 
with the allometric laws in Figure IC^ 



26 



nal flow. Therefore, before we fit tlie data clouds by using a line on NM3 and 
NM4, we actually left the data under the lines as outliers. In this way, we can 
estimate the right 7s for NM3 and NM4 which are close to 1. By comparing 
Figure IC.8I and Figure IC.6t we know that the relationship between 7 and rj 
is also suitable for these null models. 

Appendix D. Flow Adjusting Algorithm 

In the main text Section 13. 4^ we apply the so called "Flow Adjusting 
Algorithm" to perturb the original flow network. In this section, we will 
introduce the detailed steps of this method. 

The main purpose is to conform the adjusted flow network to satisfy 
two conditions (1) dissipation law with given exponent, and (2) flow balance 
requirement without changing the topology and relative weight of each flow. 
Or we can express it as a mathematical problem: to find a solution of the 
following equations system. 

N 

< ^ VzG[l,iV] (D.l) 

. i=o j=i 

Where c and 7 are given constants, fijs and DiS are variables. Therefore 
we have totally 2N equations but E + N variables. For normal flow networks, 
E > N, so we should have infinite solutions to the Equation ID.ll However, 
solving these equations are hard because they are non-linear. 

The flow adjusting algorithm is an approximate algorithm to solve these 
equations. At first, we have a nodes set O, a set of stop criterions. Initially, 
we set time step t = 1, O^^-* = {1, 2, ■ ■ -, N}, the flux matrix as the original 
flow network F^^^ = F, where the superscripts on O and F are the current 
time step. The algorithm will repeat the following steps: 

(1). For any node i in 0^^\ the algorithm needs the current out flows 
from i, {f-j^lf-j > 0}. Solve the equation for xf 

N 

i=i 
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i.e., the new total influx to i. 

(2) . Assign Xi to all incoming edges to i proportionately, set 

At) 
Z^j / ji 

and 

^(t+i) ^ 

(3) . Add all input nodes js of i into set O'-*"^^-* and delete z from O*^*^^-*; 

(4) . Set t = t + 1, repeat the previous steps until the stop criterions, 
which include the total running time being smaller than a given number, the 
dissipation law exponent 7 and for the new flux matrix F'-*-' being close 
to the wanted values, are satisfied. 

This algorithm works once the network is connected (which means there 
is at least one path from to every node i). For most of networks, the 
algorithm can converge to a network pertaining the significant dissipation 
law on the given 7 exponent. However, it may oscillate on some topologies, 
especially for random networks. Further approvement of this algorithm will 
leave for the future works. 



Appendix E. 7 and rj Relationships for Modeled Networks 

To better understand how 7 correlates with r], we will study several special 
modeled networks in this section. 

Appendix E.l. Minimum Spanning Tree 

We now consid e r a sp ecial case: minimum spanning tree introduced by 



_spe 

Frank and Murrelll ( 120051 ). By controlling two parameters 6 and /3, we can 



generate variant trees with different basal species ratio (controlled by /3) and 
maximum trophic level (controlled by 9). 

The tree's construction process is as follows. Let's consider an ecological 
community with S different species, in which a hypothetic food web (tree 
structure) will be built. At first, we select [3S species as the basal species 
at the first trophic level. And at each time, a new species j is added to the 
minimum spaning tree, j will select a node i as its unique prey according to 
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the probability: 

Pij - ^ —0, (E.l) 

where, tj is j's trophic level + 1 (i.e., the depth of j in the tree), T is the 
set of species which are already in the spanning tree, and ^ is a parameter 
to control the attachment preference of the new node on depth. If 9 is large, 
the new node may attache to the position being close to the root. 

After a tree is constructed, we will assign random values in the original 
flux matrix and then apply the "Flow Adjusting Method" on it. In this 
way, we can investigate the influence of both tree's structure {6 and (3) and 
dissipation exponent 7 on the allometric scaling exponent rj. 

From Figure IE.9t we found at first both the dissipation law and network 
structure can affect the allometric scaling. However, r] depends more on 7 
than /3 and 9 because 77 will change with 7 intensively. When 7 is given, 
we can obser ve the similar trend of depends on {3 and as the results 



introduced by iFrank and Murrelll (120051 ). 



Interestingly, when 7 is set to be 0, each node's dissipation is a constant, 
this corresponds to Garalaschelli's approach's assumption (see Figure [U^d)). 
And the exponent r] derived by our algorithm is exactly same as the result 
derived by Garlaschelli's ap proach on the same tree. And the dependence 



of 77 on /3 and 6 is same as (iFrank and Murrelll . 120051 ). So, our method can 



recover Garlaschelli's method on spanning trees. 

Appendix E.2. Random Network 

We also test the "Flow Adjusting Algorithm" on random networks. The 
results are shown in Figure IE. 101 

We only show the results of random networks with 50 nodes but differ- 
ent number of edges because the FAA is hardly to converge on the random 
networks with large number of edges. When the algorithm cannot get a final 
result after a given number of time steps (200), we have to regenerate a new 
random network with the same number of nodes and edges. It is interesting 
to observe that when the number of edges is large, the responding curve of 
?7 on 7 is very different from the ones in Figure [31 ?7S are always very small, 
where rj almost gets a peak when 7 approaches 1. Therefore, we know the 
topological structure does affect the allometric exponent. 

Furthermore, we generate random networks based on a minimum span- 
ning tree by adding a{N{N — l)/2 — A^) additional edges randomly. When 
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Figure E.9: 77 and Change with 9,/3 and 7 
For each combination of parameters, we generate 10 minimum spanning trees to 
get the average value of ij. The species number S = 100 in ah simulations 
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Figure E.IO: r\ and i?^ Change with 7 for different random networks 
All rys are averaging for 10 random networks 

a = 0, the network is an MST, while when a = 1, it is a complete graph. 
Therefore, we can observe how the corresponding curve of 77 on 7 changes 
when the network structure changes from a tree to a random network by 
tuning a. 

From Figure lE.llt the dependance of on 7 when a = 0.4 is different 
from the random networks with the same connectance but similar to the ones 
of empirical ecological networks. That is because the former random network 
is generated based on a minimum spanning tree which can be viewed as its 
backbone. Therefore, the spanning tree backbone is a key ingredient to the 
allometric exponent. Additionally, in Figure [E.llt all the random networks 
based on spanning tree with different a have a peak on 77 when 7 = 1. 

According to these experiments, we know both network structure and 
dissipation law exponent can influence allometric exponents. But the expo- 
nent 7 is more important than the structure. And the shape of backbone 
spanning trees can change the shape of the relationship between rj and 7. 

Appendix F. Flovi^ Adjustment for Empirical Ecological Flow Net- 
works 

"Flow Adjusting Algorithm" is applied to collected empirical ecological 
flow networks as shown in Figure |3] on Mondego as an example. In this 
section, we will show the results for other networks and discuss the technique 
details. 
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Figure E.ll: rj and Change with 7 and a on the Random Networks based on Spanning 
Trees 

All rjs are average results of 10 networks, the backbone trees are of the parameters 
N = 100, /3 = 0.5,6' = 5. 

The basic idea of "Flow Adjusting Algorithm" is to tune various flows on 
the network until the dissipation law exponent is close to the wanted value. 
When we apply this method on the empirical ecological networks, we adjust 
the flows until either (1) the adjusted dissipation law exponent 7' satisfies 
I7' ~ 7*1 < 0.01, where 7* is the wanted exponent and (2) the R'^ of the 
dissipation power law on the adjusted flows should satisfy R'^ < i?^, where 
R^ is the Rsquare of the dissipation law for the original flow network; or 
(3) The running time steps are larger than 500. Because the algorithm may 
diverge, whence the first requirement may not be satisfied, we have to stop 
the algorithm within a finite time steps and retrieve one of the best network 
as the output. 

In Figure IF.12t we show results for applying FAA on all collected ecologi- 
cal networks. The red diamond and green stars are the 7 and rj combinations 
for the original networks, and the red disk and green squares stand for the 
ones for FAA results on the original 7 as the designed exponent. If the 
original flow network satisfies the dissipation law and flow balance condition 
perfectly, then the adjusted r] value should be similar with the value of the 
original network (which indicates that the red disks(green squares) should 
overlap with the red diamonds (green stars)). However, we observe that the 
markers for some networks do not overlap which means the corresponding 
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Figure F.12: 7 and rj Relation Adjusted by "Flow Adjusting Algorithm" on Collected 
Ecological Networks 

Each plot shows results of "Flow Adjusting Algorithm" on two empirical networks 
(black and blue curves), the red diamond and disk stand for the 7,7/ combinations 
for the original exponents and the adjusted results by "Flow Adjusting Algorithm" 
on the original 7 for the first ecological network; and the green star and square are 
for the second ecological network. 
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Figure F.13: Dissipation Scaling Law for Original and Balanced Flow Networks of Baydry 
and Rhode 



original networks do not satisfy that two conditions perfectly. 

Let's see Baydry network as an example. Although the original flow net- 
work possess a very good dissipation law as shown in the left-top plot in 
Figure IF.13t this network is not balanced. Because balance is a basic re- 
quirement of the algorithm for allometric scaling law mentioned in Section 
we have to balance Baydry network at first by the method in Section 



Appendix A[ However, the balanced flow network always has different dis- 
sipation law (compare the left-top plot to the right-top one in Figure lF.1311 . 
Therefore, the original t] and the one adjusted by FAA are different (The 
first plot red disk and diamond in Figure [F.12p . 

Another example network which is balanced but not be of good dissi- 
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pation law is Rhode as shown in the left-bottom and right-bottom plots of 
Figure IF. 131 We observe that the dissipation scaling law is not significant 
(i?^ = 0.737) for Rhode network. Therefore, the original and adjusted expo- 
nents do not overlap (The first plot green star and square in Figure [F.12p . 
Furthermore, FAA cannot obtain a convergent result on Rhode network, 
that is the reason why the green star and square have different horizontal 
coordinates. 

Because the negative relationship between 7 and 77 is significant only if 
the flow network satisfies (1) a significant dissipation law and (2) the flows 
are balanced, we adjust the flows of the empirical networks to conform the 
dissipation law with the original dissipation law. In this way, we believe 
the noise contained in original data can be eliminated, so that more exact 
allometric exponents can be computed. The red solid curve in Figure |3] shows 
the original 7 exponents and adjusted r] exponents. 
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